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ABSTRACT 

We present radial velocities for 14 stars on the California & Carnegie Planet 
Search target list that reveal new companions. One star, HD 167665, was fit with 
a definitive Keplerian orbit leading to a minimum mass for the companion of 
50.3 Mj up at a separation from its host of ~ 5.5 AU. Incomplete or limited phase 
coverage for the remaining 13 stars prevents us from assigning to them unique 
orbital parameters. Instead, we fit their radial velocities with Keplerian orbits 
across a grid of fixed values for M sin i and period, P, and use the resulting xl 
surface to place constraints on M sin i, P, and semimajor axis, a. This technique 
allowed us to restrict M sin i below the brown dwarf - stellar mass boundary 
for an additional 4 companions (HD 150554, HD 8765, HD 72780, HD 74014). If 
the combined 5 companions are confirmed as brown dwarfs, these results would 
comprise the first major catch of such objects from our survey beyond ~ 3 AU. 

Subject headings: binaries: spectroscopic — planetary systems — stars: low- 
mass, brown dwarfs — techniques: radial velocities 



1. Introduction 

The California & Carnegie Planet Search (CCPS) has successfully identified over 120 



planetary companions to stars in the solar neighborhood with the Doppler technique (IButler et al 
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20061 ). While the focus of the survey is the search for planets, the presence of brown dwarfs 
and low-mass stellar companions in the radial velocity data serves as another avenue for 
study Brown dwarfs in particular have emerged as intriguing research targets as they oc- 
cupy the region of the mass function between stars and planets with mass boundaries of ~ 
80 Mj up and ~ 13 Mj up respectively. Unlike stars, which avoid collapse by gas pressure from 
the high temperatures and densities in their cores that ignite hydrogen fusion, brown dwarfs 
are supported from collapse by electron degeneracy pressure. For this reason, the minimum 
mass necessary to sustain hydrogen burning is taken to be the upper mass limit for brown 
dwarfs. The lower mass boundary for brown dwarfs is determined by the deuterium burning 
limit. The precise values for these transition masses are not concrete as they depend on 
parameters such as metallicity. 

With their larger imparted reflex velocity on the primary star, brown dwarfs are more 
amenable to detection than planets. However, one significant result from the first ~ 5 
years of our radial velocity survey is the lack of detected companions orbiting within ~ 
3 AU of solar-type stars with minimum mass, M s in i, in the brown dwarf regime - a 



phen omenon seen in earlier radial velocity surveys (ICampbell et al 



1993 ) and commonly referred to as the 'brown dwarf desert' ( Marcy fc Butler I2OO0I ). Planets 



1988 



Murdoch et al. 



however , are more routinely found in our survey with their mass function rising towards lower 
masses ( IMarcy et al.ll2005al ). Such a dichotomy would seem to indicate distinct formation 
mechanisms for these two classes of substellar objects. 

The lack of brown dwarf orbital companions arou nd Sun-like stars has been studied in 
detail by several groups. iGrether fc Lineweaverl (120061 ) combined data from the major radial 
velocity surveys for those stars observed over a time baseline of 5 years, resulting in uniform 
detection of all companions more massive than Jupiter within ~ 3 AU of the host star. Their 
results confirm the rising mass function for planetary companions towards lower masses and 
reveal a rising mass function for stellar companions towards higher mass, effectively leaving 
a valley around the brown dwarf desert. 

Despite the clear deficit of brown dwarf mass companions in the orbital parameter 
space currently acce ssible by the major radial velocity surveys, there are often exceptions. 
Nidever et al.l (120021 ) reported 7 companions with limited radial velocity phase coverage with 
M sin i below the hydrogen burning limit. Since that publication however, subsequent obser- 
vations leading to more complete ph ase coverage have pushed at least 1 of those companions 
into the stellar regime (HD 65430). lEndl et al.l (120041 ) also report the discovery of a brown 
dwarf candidate in the desert with minimum mass M sin i = 26 Mj up , and period P = 798 
days. However, Hipparcos astrometry measurements do not constrain the companion's mass 
to be entirely in the brown dwarf regime. 
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One explanation put forth for the presence of a desert at close separations is the merger 
of brown dwarfs with their host s tars due to orbital migrati on caused by interactions with 
an evolving protoplanetary disk (lArmitage fc Bonnelll |2002| ). The theory leaves open the 
possibility of finding br own dwarfs as companion s to solar type stars at larger separations. 
While the simulations of iMatzner fc Levin! (120051 ) also confirm the exist ence of a desert, they 
argue that brown dwarfs are not likely to form in disks as suggested by lArmitage fc Bonnell 
( 120021 ) . The argument removes a wider range of orbit al separations available to brown dwarfs 
and is strengthened by the coronagraphic survey of McCarthy fc Zuckermanl (120041 ) which 
reveals a desert that extends beyond ~ 1000 AU. 

Here we present radial velocities revealing 14 new companions from the CCPS survey. 
Some of the inferred companions are examples of potential brown dwarf candidates in the 
desert. For one star, HD 167665, we were able to fit a well-determined Keplerian orbit. 
The remaining 13 stars lack the phase coverage required to obtain a single set of orbital 
parameters. For these targets, we explore xl f° r Keplerian fits across a grid of fixed values for 
M sin i and P, thus providing constraints on these parameters. This technique is discussed 
further in § 13.31 In §|2]we discuss our observations and technique for radial velocity extraction 
from the spectra. Stellar properties and orbital parameters for new companions are presented 
in § [31 We discuss alternative means for constraining these orbits and offer concluding 
remarks in § HJ 



2. Observations 



Three telescopes and two high resolution echelle spectrometers are currently employed 
by the CCPS in the Northern H emisphere. The Keck I 10-m telescope is used in conjunction 
with HIRES (jVogt et al.lll994l ). while either the 3-m Shane or the 0.6- m Coude Au xiliary 
Telescope (CAT) at Lick Observatory feed the Hamilton spectrometer (jVogtl 119871 ). The 
spectrometers provide resolving powers of R ~ 60, 000 and 55,000 for HIRES and the Hamil- 
ton respectively. The CCPS also uses the AAT and Magellan telescopes for coverage of the 
Southern Hemisphere. 



Our radial velocity extraction technique, discussed at length in iButler et al.l (119961 ) and 
summarized below, allows us to obtain high precisio n relative radial velo cities with uncer- 
tainties of 2.5-4 m s _1 for the majority of our stars (IMarcy et al.ll2005bl ). After traversing 
the optics of the telescope, light from the star passes through a glass cell of iodine gas before 
entering the slit of the spectrometer. As a result, thousands of iodine absorption lines are 
imprinted on the observed stellar spectrum thus serving as a wavelength fiducial. A model 
stellar spectrum consisting of a high signal-to-noise template stellar spectrum multiplied by 
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the predetermined transmission function of the iodine cell is then convolved with the PSF 
of the spectrometer and fit to the observed stellar spectrum. Among the many parameters 
to the fit is the Doppler shift of the template stellar spectrum. 



The target list of stars for our survey is drawn from the Hipparcos catalogue (jESAlll997l ) 
with the cr iteria that they have B — V > 0.55. For a detailed discussion on our selection 



criteria see 



stars can be found in lWright et al.l (120041 ) 



Marcv et al.l (2005b). Select stellar properties for the majority of our program 



Observations for the 14 stars with companions presented here were obtained as far back 
as 1996. Because these companions are more massive than the planets found in our survey, 
they impart a larger radial velocity amplitude on the primary star. As a consequence, we 
typically obtain spectra for the host star with shorter integration times resulting in larger 
photon-limited velocity precision (~ 1.5-15 m s _1 ). 



Results 



3.1. Stellar Properties 



Properties for stars with companions presented in this paper are given in Table [TJ 
Columns (1) and (2) indicate the HD and Hipparcos catalogue identification numbers re- 
spectively. The spectral type in column (3) is that given by SIMBAD. When available, 
the luminosity class was also take n from SIM BAD, otherwise the derived stellar parameters 
given in the table were used with ICoxl (120001 ) to provide a best guess. All of the host stars 
are dwarfs with the exception of four subgiants (HD 211681, HD 215578, HD 29461, HD 
8765). The distance (d) and absolute visual magnitude (My), given in columns (4) and 



(5) respectivel y, are taken 



obtained from 



Takeda et a 



tained from Girardi et al. 



r om H ipparcos. Masses for the host stars in column (6) were 
(120071) with the exception of HD 215578, whose mass was ob- 
J2002h . The next four columns indicate surface temperature (7), 



surface gravity (8 ) , met allicity (9), and rotation velocity (10), which were obtained from 
Valenti &: Fischerl (120051 ). The primary observing facility is given in column (11). 



3.2. HD 167665 

Orbital parameters for HD 167665 are given in Table [2j These were derived from times 
of observations, radial velocities and instrumental uncertainties that are given in Table [3j 
Radial velocities with best-fit Keplerian orbit are shown in Figure [TJ Instrumental uncer- 



tainties on the individual radial velocity measurements are also plotted in this figure but are 
in most cases too small to see at the given scale. A Levenberg-Marquardt algorithm was 
implemented to determine the best Kepl erian fit to th e observ ed data. Uncertainties for the 



orbital parameters were determined as in lMarcy et al.l (l2005bl ). The residuals from the best- 



fit Keplerian orbit were treated as samples of the uncertainty due to the combined Doppler 
error and photospheric jitter. One hundred mock velocity sets were then created where a 
new mock velocity was constructed by adding the original observed velocity to a randomly- 
picked residual. The final uncertainties for the orbital parameters as reported in Table [2] 
are given by the standard deviation of the best-fit parameters to the 100 mock velocity sets. 
The uncertainty for M sin i was calculated by propagating the errors for the relevant orbital 
parameters through the mass function equation. Note that the quoted uncertainty for M sin 
i does not reflect the uncertainty associated with the host mass, which is typically ~ 5-10%. 
The error for the minimum semimajor axis, a m ; n , was computed by propagating the relevant 
errors through Kepler's 3rd Law. 

The companion to HD 167665 orbits with a period of ~ 12 yr at a separation from its 
host of ~ 5.5 AU. It has a minimum mass firmly established in the brown dwarf regime, with 
M sin i = 50.3 Mj up . The true mass depends on the inclination of the orbit with higher 
inclinations leading to larger true masses. If we assume the orbit is randomly oriented on 
the unit sphere, there is a 79% chance that the inclination favors an object in the brown 
dwarf mass range. We take 82 Mj „n as the mass bounda ry between stars and brown dwarfs 



by interpolating between tracks in iBurrows et al.l (119971 ). 



3.3. Stars with Limited Phase Coverage 

Thirteen of the 14 stars presented in this paper were not observed long enough to realize 
a full orbit. Any attempt to present a single set of orbital parameters for them would be 
misleading as any number of a large family of solutions could result in an acceptable fit. 
However, if the radial velocity data shows enough curvature it may be possible to put limits 



on th e period (P) and minimum mass (M sin i). Following the analysis in IWright et al. 



(120071 ). we examine the reduced chi-square statistic, xti resulting from Keplerian fits with 
fixed values of M sin i and P. Typically the semiamplitude, K, is the parameter of interest 
when detecting spectroscopic companions. In order to place constraints on the minimum 
mass we transform from (P, e, K) to (P, e, M sin i) coordinates through the mass function 
equation, f(M): 

. M 3 sin 3 i PK 3 (1 - e 2 )l 
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For each of our limited phase coverage targets, we create a grid of values for M sin i 
and P, both in log-step. These are the fixed values with which the Keplerian fitting routine 
is executed while allowing e, u>, T P and 7 to vary. Again, we use a Levenberg-Marquardt 
algorithm to arrive at a solution, which is in almost all cases the global minimum given the 
fewer free parameters to the fit and because of the velocity precision at these large amplitudes. 
A smoothly varying xl surface on the M sin i - P grid provides further evidence that we 
are tracing the global minimum. 

An initial grid with a large range in period and minimum mass was constructed to 
determine the region of M sin i - P space to iterate on. In Figures 121 - [T41 we show the more 
focused grids of xl m M sin i — P space, with contour levels of 1, 4, and 9 from the minimum 
(xlmm) corresponding to 1, 2, and 3-sigma confidence levels. One hundred logarithmic steps 
in M sin i and P are used for each focused grid with step sizes dependent on the range 
displayed on the axes in Figures 2-14. Contours of eccentricity are shown in these figures in 
increments of 0.2. Also presented in the figures are the radial velocity measurements for each 
star along with the Keplerian orbit corresponding to the best-fit on the M sin i - P grid 
(indicated by a cross). The best-fit xl quoted in these figures includes the contribution from 
the estimated photospheric jitter, however we caution the reader to note the uncertainty in 
this value due to the uncertain nature of the jitter. The M sin i - P grids were run without 
the jitters for this reason and because the results are insensitive to the jitter levels assumed 
at these large amplitudes. As a test, we find that when including jitter, the minimum xl 
in the M sin i- P grid corresponds to the best-fit solution we obtain when excluding jitter. 
T he jitter s given in the figures were derived from recent improvements to the jitter algorithm 



of IWrightl (120051 ). We note that while it is often customary to remove linear velocity trends 
before performing the Keplerian fits, we refrain from doing so here because the limited phase 
coverage prevents a well-determined linear solution. 

Using the xl — xlmm + 4 contour as a constraint, we can place limits on the period and 
minimum mass. We also implement an additional eccentricity constraint by noting that for 
those s pectroscopic binaries w ith P > 10 yr in the Ninth Catalogue of Spectroscopic Binary 



Orbits (IPourbaix et al.ll2004j ). approximately 90% have an eccentricity e < 0.8. Therefore, 
we will assume the companions reported here have e < 0.8 unless their eccentricity explicitly 
exceeds this limit. We report the possible range in minimum mass and period along with the 
semimajor axis, a, in Table HI Note that upper limits are absent when the xl — xlmm + 4 
contour extends off of the grid and the eccentricity constraint is unable to provide additional 
limits. 
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3.4. New Companions 

The orbit for the companion to HD 142229 (Figure [2]) is not well constrained. The 
location of the best-fit on the M sin % - P grid is on its upper edge (arrow), indicating that 
higher masses would lead to acceptable fits. However, a higher mass is unlikely given that 
a companion of mass greater than 1000 M Jup (~ 1 M ) would probably have been detected 
either visually or spectroscopically. 

While the period for the companion to HD 150554 (Figure [3]) is uncertain, the minimum 
mass is well constrained due to phase coverage near the turning points. If we assume e < 0.8 
then the minimum mass for the companion is bound to the brown dwarf mass range with a 
period between 8 and 34 yr. 

HD 211681 (Figure H]) is another example of how instituting an eccentricity limit of 
e < 0.8 can help constrain the period. The companion has a period of approximately 
between 10 and 100 yr. Though the upper limit of M sin i is 102 Mj up , the lower limit does 
not exclude the companion from the brown dwarf mass regime. 

Only lower limits to M sin % and P were obtained for the companion to HD 215578 
(Figure [5]) because the xl — xl mm + 4 contour extended off of the grid. With a lower limit 
for M sin % of ~ 0.5 M , this companion is one the largest of those reported in this paper. 

Lower limits to M sin i and P are also reported for HD 217165 (Figure [6]) for the same 
reason as HD 215578. Given the curvature present in the radial velocity data, a wide range 
of minimum masses yielded acceptable fits. While the best-fit is ~ 200 Mj up , the lower limit 
for M sin i extends into the brown dwarf mass range. Note that the break in the \ 2 V surface 
between 10 and 20 yr is an artifact of proceeding to the next pixel in the grid using the 
previous pixel's solution. Traversing the grid in this fashion allows us to minimize random 
departures from the best Keplerian solution. 

The most recent radial velocity observation for HD 29461 (Figure [7j) was able to put 
strong constraints on the minimum mass and period for the companion. Assuming e < 0.8, 
the object is most likely a low-mass star with a lower limit for M sin i of ~ 90 Mj up . The 
period is between ~ 10 and 25 yr. 

The eccentricity criteria of e < 0.8 was not implemented for HD 31412 (Figure [8]) given 
its highly eccentric orbit. The companion has a minimum mass between ~ 370 and 415 
Mj up , and period between ~ 75 and 200 yr. 

The companion to HD 5470 (Figure [9]) is most likely a low-mass star. Assuming e < 0.8, 
the lower limit for M sin i is ~ 160 Mj up and the upper limit ~ 400 Mj up . The period is as 
short as ~ 13 yr and as long as ~ 170 yr, putting the companion as far away as 35 AU from 
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its host star. 

The xl — Xymin + ^ contour limits the companion for HD 8765 (Figure [TO]) to the 
brown dwarf mass regime with an upper limit for M sin i of ~ 76 Mj up . The period ranges 
between ~ 6 and ~ 18 yr. We note that HD 8765 is given an acceleration solution (G) in the 
Hipparcos catalogue indicating that a quadratic or cubic astrometric s olution was obtained 



using the position, proper motion, and parallax (ILindegren et al.lll997l ). 



The minimum mass for the companion to HD 199598 (Figure [TT]) is well-determined 
due to radial velocity measurements around the turning points, providing strong constraints 
on the semiamplitude. Assuming e < 0.8, the companion is a low-mass star with minimum 
mass between 105 and 120 Mj up and period between 25 and 85 yr. 

The companion to HD 72780 (Figure [121 is another potential brown dwarf candidate 
with a minimum mass constrained between ~ 52 and 60 Mj up , assuming e < 0.8. The period 
is less constrained, ranging between ~ 12 and 42 yr. 

The low-mass star that orbits HD 73668 (Figure [TBI is in a fairly eccentric orbit as seen 
in the radial velocity measurements. Assuming e < 0.8, the minimum mass and period are 
very well constrained. The minimum mass ranges between ~ 160 and 180 Mj up , while the 
period is bound between ~ 17 and 24 yr. We note that HD 73668 is given a component 
solution (C) in the Hipparcos catalogue indicating that a linear astrometric solution was 
obtained for it and another nearby visible component. 

Phase coverage near the turning points of HD 74014 (Figure IT4"|) puts strong limits on 
its semiamplitude and in turn the minimum mass of the companion. Assuming e < 0.8, the 
the minimum mass is between ~ 47 and 57 Mj up . The period is less constrained, with a 
lower limit of ~ 12 yr and an upper limit of ~ 77 yr. 



4. Discussion 

Of the 14 new companions to stars in the CCPS reported here, 7 have M sin i extending 
into the brown dwarf mass regime. Of those 7, 5 have upper limits for M sin i below the 
brown dwarf - stellar mass boundary (HD 150554, HD 167665, HD 8765, HD 72780, HD 
74014), making them potential brown dwarf candidates. The sin i degeneracy ultimately 
determines the true mass of these companions. 

Hipparcos data can serve as a useful tool in constraining orbital inclination and thus the 
companion mass. There are two ways in which one typically goes about doing this. The first 
method makes use of flags in the Hipparcos Double and Multiple Systems Annex indicating 
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an astrometric perturbation. One then assumes a lower limit to the astrometric perturbation 
detectable by Hipparcos, t hus providing a lower limit for the inclination (see Equation 1 of 
Pourbaix fc Arenoul (120011 )) and an upper limit for the mass of the companion. This simple 
technique should serve only as a rough guide for estimating upper limits to masses given 
that there is no consensus on what the detectability threshold is for Hipparcos. The second 
technique uses the Hipparcos Intermediate Astrometric Data al ong with spectroscopic data t o 
perform a detailed 12-parameter fit to estimate the inclination (IReffert fc Quirrenbachl|2006l ). 
Both techniques only yield meaningful results when attempting to constrain inclinations for 
those orbits with periods less than or comparable to the Hipparcos mission lifetime (~ 4 
yr). Therefore, given the long period nature of the companions presented in this paper, we 
refrain from trying to present constraints on inclination with Hipparcos data. 

Coronographic studies with adaptive optics (AO) can aid in placing further constraints 
on the masses of the brown dwarf candidates. Limiting factors in such observations include 
flux contrast ratios and the separation on the sky of the brown dwarf and the host star. 

Com bining upper l i mit M sin i estimates for our brown dwarf candidates with stellar 
ages from iTakeda et al.l (120071) , we can estimate near -IR apparent magnitudes using brown 
dwarf evolutionary models (IBurrows et al.lll997l . 120061 ). We assume the brown dwarf and the 
host star are coeval and at the same distance. The host stars have K ~ 5 - 7. Implementing 



filter transmission curves from 2MASS (jCohen et al.ll2003l ) we find the brown dwarf candi- 
dates have K ~ 17 - 21, resulting in contrasts of ~ 11 - 14 mag (~ 10~ 4 ' 6 — 10 -5 7 ). If we 
assume all of the bro wn dwarf candidates have a mass of 73.3 Mj up , the highest mass avail- 
able in the models of IBurrows et al.l (119971 ) for which log g and T e g- are tabulated, we obtain 
contrasts of ~ 10 -3 ' 9 — 10" 5 ' 3 in X-band. These contrast levels would make detectability 
with AO challenging, especially if the orbital periods and separations reside towards the 
lower limits derived from the M sin i - P grid. However, null detections would be helpful 
in providing upper limits on the masses of the brown dwarf candidates. 

In general, the relatively older ages of stars in the CCPS implies that potential brown 
dwarf companions will be quite faint. One can conclude then that imaging brown dwarfs in 
our survey within ~ 10 AU will prove to be a difficult task due to large flux contrast ratios 
at short projected separations. 

We have reported radial velocities for 14 new companions to stars on the CCPS target 
list. For HD 167665 we were able to fit a Keplerian orbit with well-determined orbital 
parameters. For the remaining stars, we explored Keplerian fits across a grid of values for 
M sin i and P. The technique allowed us to place constraints on the minimum masses, 
periods, and semimajor axes for these limited phase coverage targets. We find that 5 of the 
companions reported in this work have M sin i firmly established in the brown dwarf mass 
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range at orbital distances of ~ 4 - 18 AU. Imaging with adaptive optics could further aid 
in constraining these orbits. Confirmation of these companions as brown dwarfs would add 
more examples of them to the rare sample that exists in the brown dwarf desert. 
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Fig. 1. — Velocity vs. time for HD 167665 (dots). The solid line is the best-fit Keplerian 
orbit with P = 12.0 yr, K = 610 m s _1 , e = 0.34, and M sin z = 50.3 M Jup . 
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Fig. 2. — Left: Contours of xl i n M sin i - P space for HD 142229. The solid contour lines 
represent increases from the best-fit {xlmm) °f 1) 4, and 9. The location of the best-fit on 
the M sin % - P grid is indicated by the arrow on the upper edge at 1000 Mj up and 35.3 yr, 
and its corresponding Keplerian orbit is plotted on the right. Contours of eccentricity are 
shown as dashed lines. The xl = xl mm + 4 contour limits the minimum mass and period to: 
M sin i (M Jup ) > 153, P (yr) > 16.4. Right Velocity vs. time for HD 142229 (dots). The 
solid line is the Keplerian orbit for the best-fit (x^min) from the M sin % - P grid. 
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Fig. 3. — Left Contours of xt i n M sin i — P space for HD 150554. The solid contour lines 
represent increases from the best-fit (xl,mm) °f 1j 4, and 9. The location of the best-fit on 
the M sin i — P grid is indicated with a cross at 68.8 Mj up and 11.8 yr, and its corresponding 
Keplerian orbit is plotted on the right. Contours of eccentricity are shown as dashed lines. 
The xt = Xumin + ^ contour and e < 0.8 constraint limits the minimum mass and period to: 
65.7 < M sin i (M Jup ) < 75.2, 8.0 < P (yr) < 33.7. Right Velocity vs. time for HD 150554 
(dots). The solid line is the Keplerian orbit for the best-fit (Xumin) fr° m the M sin % - P 
grid. 
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Fig. 4. — Left Contours of xl i n M sin i — P space for HD 211681. The solid contour lines 
represent increases from the best-fit (xl,mm) °f 1j 4, and 9. The location of the best-fit on 
the M sin i - P grid is indicated with a cross at 87.8 Mj up and 27.3 yr, and its corresponding 
Keplerian orbit is plotted on the right. Contours of eccentricity are shown as dashed lines. 
The xl = xlmin + ^ contour and e < 0.8 constraint limits the minimum mass and period to: 
71.7 < M sin % (M Jup ) < 102, 10.4 < P (yr) < 106. Right Velocity vs. time for HD 211681 
(dots). The solid line is the Keplerian orbit for the best-fit (Xumin) fr° m the M sin % - P 
grid. 
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Fig. 5. — Left: Contours of xt i n M sin i - P space for HD 215578. The solid contour lines 
represent increases from the best-fit (xlmm) °f 1j 4, and 9. The location of the best-fit on 
the M sin i - P grid is indicated with a cross at 758 Mj up and 56.4 yr, and its corresponding 
Keplerian orbit is plotted on the right. Contours of eccentricity are shown as dashed lines. 
The xl — + 4 contour limits the minimum mass and period to: M sin i (M Jup ) > 523, 

P (yr) > 26.8. Right: Velocity vs. time for HD 215578 (dots). The solid line is the Keplerian 
orbit for the best-fit (xlmm) from the M sin i - P grid. 
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Fig. 6. — Left: Contours of xt i n M sin i - P space for HD 217165. The solid contour lines 
represent increases from the best-fit (xlmm) °f 1j 4, and 9. The location of the best-fit on 
the M sin i - P grid is indicated with a cross at 211 Mj up and 33.7 yr, and its corresponding 
Keplerian orbit is plotted on the right. Contours of eccentricity are shown as dashed lines. 
The xl — mm + 4 contour limits the minimum mass and period to: M sin i (M Jup ) > 45.9, 
P (yr) > 11.1. Right: Velocity vs. time for HD 217165 (dots). The solid line is the Keplerian 
orbit for the best-fit (xl mm) from the M sin % - P grid. 
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Fig. 7. — Left: Contours of xl i n M sin i - P space for HD 29461. The solid contour lines 
represent increases from the best-fit (xtmin) °f ^ 4, and 9. The location of the best-fit on 
the M sin i — P grid is indicated with a cross at 193 Mj up and 18.0 yr, and its corresponding 
Keplerian orbit is plotted on the right. Contours of eccentricity are shown as dashed lines. 



The x 



A /. 



+ 4 contour and e < 0.8 constraint limits the minimum mass and period to: 



91.2 < M sin % (M Jup ) < 224, 10.2 < P (yr) < 25.9. Right Velocity vs. time for HD 29461 
(dots). The solid line is the Keplerian orbit for the best-fit (Xumin) fr° m the M sin % - P 
grid. 
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Fig. 8. — Left Contours of xt i n M sin i - P space for HD 31412. The solid contour lines 
represent increases from the best-fit (x^ min ) of 1, 4, and 9. The location of the best-fit on 
the M sin i — P grid is indicated with a cross at 395 Mj up and 127 yr, and its corresponding 
Keplerian orbit is plotted on the right. The xl — xlmm + 4 contour limits the minimum 
mass and period to: 368 < M sin % (M Jup ) < 415, 75.7 < P (yr) < 205. Right Velocity vs. 
time for HD 31412 (dots). The solid line is the Keplerian orbit for the best-fit (x^mm) fro m 
the M sin i — P grid. 
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Fig. 9. — Left: Contours of xl i n M sin % - P space for HD 5470. The solid contour lines 
represent increases from the best-fit (xtmm) °f ^ 4, and 9. The location of the best-fit on 
the M sin i - P grid is indicated with a cross at 243 Mj up and 25.5 yr, and its corresponding 
Keplerian orbit is plotted on the right. Contours of eccentricity are shown as dashed lines. 
The xl = xlmin + 4 contour and e < 0.8 constraint limits the minimum mass and period to: 
163 < M sin i (M Jup ) < 425, 13.0 < P (yr) < 167. Right: Velocity vs. time for HD 5470 
(dots). The solid line is the Keplerian orbit for the best-fit (Xumin) fr° m the M sin % - P 
grid. 



-22 - 




r " n Time(yr) 

Fig. 10. — Left: Contours of xl i n M sin i — P space for HD 8765. The solid contour lines 
represent increases from the best-fit (xl,mm) °f 1j 4, and 9. The location of the best-fit on 
the M sin i — P grid is indicated with a cross at 51.5 Mj up and 7.69 yr, and its corresponding 
Keplerian orbit is plotted on the right. Contours of eccentricity are shown as dashed lines. 
The xl = xlmin + ^ contour and e < 0.8 constraint limits the minimum mass and period to: 
43.0 < M sin i (M Jup ) < 75.7, 6.3 < P (yr) < 18.1. Right Velocity vs. time for HD 8765 
(dots). The solid line is the Keplerian orbit for the best-fit (xtmin) fr° m the M sin % - P 
grid. 
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Fig. 11. — Left: Contours of xl m M sin i — P space for HD 199598. The solid contour lines 
represent increases from the best-fit (xl,mm) °f 1) 4, and 9. The location of the best-fit on 
the M sin i - P grid is indicated with a cross at 112 Mj up and 59.0 yr, and its corresponding 
Keplerian orbit is plotted on the right. Contours of eccentricity are shown as dashed lines. 



The x 



+ 4 contour and e < 0.8 constraint limits the minimum mass and period to: 



105 < M sin i (M Jup ) < 120, 25.6 < P (yr) < 85.3. Right Velocity vs. time for HD 199598 
(dots). The solid line is the Keplerian orbit for the best-fit (Xumin) fr° m the M sin % - P 
grid. 
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Fig. 12. — Left: Contours of xl i n M sin i — P space for HD 72780. The solid contour lines 
represent increases from the best-fit (xl,mm) °f 1j 4, and 9. The location of the best-fit on 
the M sin i — P grid is indicated with a cross at 55.3 Mj up and 17.5 yr, and its corresponding 
Keplerian orbit is plotted on the right. Contours of eccentricity are shown as dashed lines. 
The xl = xlmin + ^ contour and e < 0.8 constraint limits the minimum mass and period to: 
52.4 < M sin i (M Jup ) < 60.4, 11.8 < P (yr) < 41.7. Right Velocity vs. time for HD 72780 
(dots). The solid line is the Keplerian orbit for the best-fit (^ min ) from the M sin % - P 
grid. 
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Fig. 13. — Left: Contours of xl i n M sin i — P space for HD 73668. The solid contour lines 
represent increases from the best-fit (xl,mm) °f 1) 4, and 9. The location of the best-fit on 
the M sin i - P grid is indicated with a cross at 164 Mj up and 20.4 yr, and its corresponding 
Keplerian orbit is plotted on the right. Contours of eccentricity are shown as dashed lines. 



The x 



A /, 



+ 4 contour and e < 0.8 constraint limits the minimum mass and period to: 



160 < M sin % (M Jup ) < 181, 17.1 < P (yr) < 24.1. Right Velocity vs. time for HD 73668 
(dots). The solid line is the Keplerian orbit for the best-fit (Xumin) fr° m the M sin % - P 
grid. 
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Fig. 14. — Left: Contours of xl m M sin i — P space for HD 74014. The solid contour lines 
represent increases from the best-fit (xl,mm) °f 1j 4, and 9. The location of the best-fit on 
the M sin i - P grid is indicated with a cross at 49.7 Mj up and 19.3 yr, and its corresponding 
Keplerian orbit is plotted on the right. Contours of eccentricity are shown as dashed lines. 
The xl = Xumin + ^ contour and e < 0.8 constraint limits the minimum mass and period to: 
46.6 < M sin i (M Jup ) < 57.3, 11.9 < P (yr) < 77.4. Right Velocity vs. time for HD 74014 
(dots). The solid line is the Keplerian orbit for the best-fit (^ min ) from the M sin % - P 
grid. 
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Table 1. Stellar Properties 









d 


M v 










V sin i 




HD 


HIP 


Spectral Type 


(pc) 


(mag) 


(Mo) 


(K) 


log g 


[Fe/H] 


(km s -1 ) 


Observatory 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


142229 


77810 


G5V 


40.5 


5.04 


1.06 


5930 


4.50 


0.050 


5.5 


Keck 


150554 


81662 


F8V 


44.8 


4.43 


1.07 


6010 


4.29 


0.0 


2.7 


Keck 


167665 


89620 


GOV 


29.7 


4.04 


1.09 


6115 


4.22 


-0.17 


5.2 


Keck 


211681 


109169 


G5 IV/V 


70.9 


3.85 


1.31 


5839 


4.32 


0.45 


2.6 


Keck 


215578 




K0 IV/V 


43.2 


4.33 


1.04 


5136 


3.37 


0.43 


3.9 


Keck 


217165 


113438 


GOV 


43.7 


4.47 


1.05 


5936 


4.39 


0.010 


3.2 


Keck 


29461 


21654 


G5 IV/V 


48.1 


4.53 


1.10 


5913 


4.52 


0.25 


0.9 


Keck 


31412 


22919 


F8V 


36.0 


4.24 


1.13 


6096 


4.37 


0.050 


2.6 


Keck 


5470 


4423 


GOV 


67.8 


4.19 


1.13 


5966 


4.32 


0.24 


1.4 


Keck 


8765 


6712 


G5 IV/V 


75.1 


3.77 


1.04 


5590 


4.23 


0.19 


2.7 


Keck 


199598 


103455 


GOV 


33.2 


4.32 


1.11 


6022 


4.40 


0.070 


1.9 


Lick 


72780 


42112 


F8V 


52.5 


3.87 


1.23 


6266 


4.37 


0.15 


7.0 


Lick 


73668 


42488 


G1V 


35.7 


4.50 


1.05 


5941 


4.41 


0.0 


3.2 


Lick 


74014 


42634 


KOV 


33.5 


4.96 


0.968 


5605 


4.43 


0.26 


2.4 


Lick 
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Table 2. Orbital Parameters for HD 167665 



Parameter HD 167665 

P (days) 4385 ± 64 

T p (JD-2,440,000) 8140 ± 70 

e 0.337 ± 0.005 

u> (deg) 225.0 ± 1.1 

K (m s- 1 ) 609.8 ± 3.1 

amin (AU) 5.47 ± 0.05 a 

M sin i (Mj up ) 50.3 ± 0.4 a 

rms (m s -1 ) 6.45 

N obs 23 



a Uncertainty does not include con- 
tribution from host mass. 



Table 3. Radial Velocities for HD 167665 





RV 


Uncertainty 


JD-2,440,000 


(m s- 1 ) 


(ms- 1 ) 


10283.918 


436.87 


2.10 


10604.041 


329.14 


1.67 


10665.871 


295.83 


1.84 


10955.098 


170.89 


2.29 


11013.928 


151.26 


2.50 


11069.828 


101.16 


2.10 


11313.053 


-33.63 


2.18 


11341.895 


-41.21 


1.95 


11367.816 


-61.46 


2.08 


11410.805 


-100.66 


2.57 


11439.715 


-111.51 


1.81 


11679.030 


-281.80 


1.86 


11702.962 


-282.59 


2.23 


12007.088 


-522.40 


2.00 


12095.949 


-588.79 


2.42 


12162.750 


-599.04 


2.09 


12390.040 


-569.74 


2.51 


12488.861 


-482.93 


2.36 


12574.701 


-373.79 


2.19 


13072.152 


347.96 


2.20 


13238.776 


476.40 


1.93 


13604.757 


590.64 


2.75 


13838.129 


599.04 


2.76 
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Table 4. Constraints for Limited Phase Coverage Targets 



M sin i P a 

Wjup) (yr) (AU) 



Star min max min max min max 



(1) (2) (3) (4) (5) (6) (7) 



HD 142229 


153 




16.4 




6.9 




HD 150554 


65.7 


75.2 


8.0 


33.7 


4.2 


11 


HD 211681 


71.7 


102 


10.4 


106 


5.3 


25 


HD 215578 


523 




26.8 




10 




HD 217165 


45.9 




11.1 




5.1 




HD 29461 


91.2 


224 


10.2 


25.9 


5.0 


9.6 


HD 31412 


368 


415 


75.7 


205 


20 


40 


HD 5470 


163 


425 


13.0 


167 


6.0 


35 


HD 8765 


43.0 


75.7 


6.3 


18.1 


3.5 


7.1 


HD 199598 


105 


120 


25.6 


85.3 


9.3 


21 


HD 72780 


52.4 


60.4 


11.8 


41.7 


5.6 


13 


HD 73668 


160 


181 


17.1 


24.1 


7.1 


8.9 


HD 74014 


46.6 


57.3 


11.9 


77.4 


5.2 


18 



